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Abstract. The problem of reconstructing the drift of a diffusion in R'', d > 2, from 
the transition probability density observed outside a domain is considered. The solu- 
tion of this problem also solves a new inverse problem for a class of parabolic partial 
differential equations. This work considerably extends W in terms of generality, both 
concerning assumptions on the drift coefficient, and allowing for non-constant diffu- 
sion coefficient. Sufficient conditions for solvability of this type of inverse problem for 
d = 1 are also given. 



1. Introduction 
Let (xj,P^) be the weak solution of the stochastic equation 

(1) Xt = X + c{xs)ds+ / y/a{xs) dwg 

Jo Jo 

where c is a measurable vector field on M'^, d > 2, a is symmetric positive-definite ma- 
trix, and w. is a standard d-dimensional Brownian motion. Denoting with p{x, t, y) the 
transition probability density of (xt,P^), our aim is to reconstruct c from the observa- 
tion of p{x, t, y) for (x, y) G A'^ x A'^, t G [0, T], where A is a bounded domain in and 
T > 0. We point out that this is an entirely different problem from the ones which are 
usually studied in filtering theory (or statistics of processes where the observations are 
partial with respect to time instead of space) . 

The problem we handle in the present paper was introduced and solved in 2^ under 
the assumptions that o(x) is constant, c = V99 with Lp G C'^iW^), and under growth con- 
ditions on c such that a strong solution to exists with infinite lifetime. The results 
were rediscussed in 0, where more detailed proofs can be found, and an approach to 
the case of diagonal diffusion through a random time change is sketched. We improve all 
these results in several directions. In particular, we impose only integrability assump- 
tions on if, we require only the existence of weak solutions instead of strong ones, and 
we allow A to be unbounded. Moreover, we consider much more general classes of SDKs 
with variable diffusion coefficients. We also obtain a partial solution of the problem in 
the case of the drift not being a gradient field and for one-dimensional diffusions. In 
their full generality, however, these latter problems were and remain unsolved (to the 
best of our knowledge). 

The problem at hand admits a purely analytic interpretation. Namely, Kolmogorov's 
classical work implies that p{x, t, •) solves the parabolic partial differential equation 
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(PDE) 

(2) (1 = ^-^' 

where L* ^ is the formal adjoint of the operator 

Lc,au{x) = ^ tr {a{x)D'^u{x)) + {c{x),Du{x)) 

(here and throughout the paper we adopt the convention of summation over repeated 
indices). Therefore our problem can be reformulated as an inverse problem for PDEs: 
from the knowledge of the solution of ((2)) in A'^ alone, determine the coefficient c of the 
first-order term (as a function on M.'^). To the best of our knowledge the problem has 
not been addressed in the literature on inverse problems for PDEs (see e.g. (Hj, |15j). 
and our solution can also be seen as a probabilistic solution to this analytic problem. 

The paper is organized as follows: we collect basic assumptions, definitions, and some 
known results in section |2l In section 3 we derive a representation formula for the 
transition probabilities of diffusions whose generators satisfy certain conditions (this 
class, in particular, contains distorted Brownian motion - see e.g. jHI). We also prove 
some consequences of this representation formula that allow us to reconstruct from the 
transition densities of a diffusion a function of its drift coefficient (in section 4), and 
eventually its drift (in section 5). Section 6 deals with extensions such as unbounded A 
and one-dimensional diffusions. 

2. Preliminaries 

Unless otherwise stated, we shall work under the following standing assumptions: 

(i) the operator Lc^a is uniformly elliptic, i.e. 

{a{x)tO>m^ Vx,eGM^ 

for some constant 6 > 0; 

(ii) equation admits a unique weak solution in W^; 

(iii) the transition probability measures of all considered diffusions admit a continu- 
ous density with respect to Lebesgue measure. 

Remark 1. By a classical result of Zvonkin and Krylov [22]) a weak solution of equation 
exists and is unique if the coefficients are bounded and a is continuous. A detailed 
study of conditions implying that a (generalized) diffusion has continuous probability 
densities can be found in TO^ (see also jlH] and [HI). 



We shall make use of some Banach spaces. We denote by := LP(IR , M), for 
p G [1, -|-oo[, the space of functions / : M'^ ^ M such that 



\f\l.--= I \f{x)\Pdx<^. 



For s G M, let RP'' = (1 - /\)-'/^LP be the usual space of Bessel potentials on M . 
The norm in H^'^ can be taken to be equivalent to the one of W^'"^, the usual Sobolev 
space of function with (generalized) derivatives up to order 2 in L^. Hf^^ is the space of 
functions / such that f( G H^'^ for all C € Cq°{W^). The following Sobolev embedding 
theorem holds: / G Hf^^(U.'^) implies that / G C"(M'^), a := s — d/p. In particular 
/ G Hf'l implies / G 
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Let us briefly recall some results from pinned diffusions connected to representation 
of transition probability densities of diffusions (see (17) and jl9j for more details). Let 
(xt,P^) be a diffusion process on M.'^, endowed with its natural filtration {J^t)t>o- 
a probability measure on a{Tt : t < T) by 

p{xt,T -t,y) 



). Define 



dF^y 



dl 



yt < T. 



The diffusion process {{xt)t<T-,Fj^^) is called a pinned diffusion, or a diffusion condi- 
tioned on xo = X and xt = y- 

The following Cameron-Martin formula for pinned diffusions was proved in |19j . 
Lemma 2. Let (xt,F^) be a diffusion process with generator 



Let c 



Lb,af = -iv{a{x)D^f{x)) + {h{x),Df{x)). 
he measurable and such that 



exp 



T 



\a '^^'^{xs)c{xs)\'^ ds 



< oo. 



Let ph and Pb+c be the transition density functions corresponding to the diffusions with 
generators L^^a ctnd Lb+c,a, respectively. Suppose that 



y 



exp 



( j\a-^'\xt)c{xt),dwt) \a-^'\xt)c{xt)\^ dt 



is continuous. Then one has 
Pb+c{x,T, y) 



(3) 



Pbix,T,y) 



exp 



(a ^^'^{xt)c{xt),dwt 



-'/^{xt)c{xt)\^dt 



3. Transition densities for a class of diffusions 



From now on we assume that the hypotheses of Lemma |21 are satisfied and that there 
exists a function tp : F., ip G H'^''^, such that Vip{x) = a^^{x)c{x). This assump- 

tion, even though it looks quite restrictive, contains many important examples. Let us 
consider, for instance, the case of unit diffusion coefficient, i.e. Markov processes with 
generator Lu = ^Au + {c, Vn). Assume that c is regular and L admits an infinitesimally 
invariant measure u{dx) = p{x) dx, i.e. that 

Lfu{dx)=0 V/gCo°°(M'^) 

(see e.g. jl] for more details on these notions). It is well known that the diffusion 



is reversible if and only if L is symmetric in L 



One can also prove that L is 



symmetric in L^(]R'',zv) if and only if 2c = V(logyo). Therefore, our problem of drift 
reconstruction can be solved for a large class of reversible diffusions. However, the 
general case of c not being a gradient field is unfortunately not within the reach of our 
method, except for some special situations discussed in the last section. 

In this section we shall specialize Lemma I^J to the class of diffusions just introduced, 
and deduce some important corollaries. 



Proposition 3. One has 

Pb+c{x, t, y) = pb{x, t, y) exp (^/;(y) - ^'(a;)) E^'^ 
where V{x) := Lb^aip{x) + ^|a^/^(x)V^/^(x)|^. 



exp 



V{xs) ds 
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Proof. By an application of Ito's lemma (as in p[B]) we get 

t r-t 

{a^^/'^{xs)c{xs),dws) + / Lb,a'4'{Xs)ds. 



Therefore 

/ {a^^^'^{xs)c{xs),dws) I \a^^^'^{xs)c{xs)\'^ ds 



il^ixt) - %Ij{xo) - I Lb^a'4^{xs)ds - ^ I \a ^/'^{xs)c{xs)\'^ ds 







il^ixt) - ip{xo) - I Lb^aip{xs)ds - ^ I |a^/^(xs)VV'(xs)p(is 



= tpixt) - tp{xo) - V{xs)ds, 
Jo 

By a simple rewriting of © we get the desired result. □ 

If 6 = 0, a = /, i.e. xt is Brownian motion, we recover a formula already obtained in 
[2], although under more regularity assumptions: 

Corollary 4. Assume that c satisfies the hypotheses of Lemma\^with 6 = and a = I. 

Setting V{x) = i(|VV'(x)p + AV'(x)), one has 

where pQ{x,t,y) is the transition probability density of Brownian bridge. 

Given a domain A C M'^ and x, y £ dA, we shall denote by (7x^)sg[o,t] the (straight) 
line joining x with y in time t, i.e. the function 

j-y = x+{y-x)^, 0<s<t. 

The following proposition will be crucial for the solution of our problem. First, let us 
observe that the elementary relation 

V{j^y) ds = t [ V{x + {y- x)s) ds 



holds. 

Proposition 5. Assume that V is bounded from below and satisfies the following prop- 
erty: 

where 7„ are continuous curves in with endpoints x, y. Then V o-j^y E L"'^(IR;M) for 
almost all x, y (with respect to Lebesgue measure) and one has 

(4) liin Pb+c{x,t,y) ^ ^^{y)-^(x) 

t-o pb{x,t,y) 

and 

(5) hm i ( log - my) - ^(x))) = - Vix + (y - x)s) ds. 
t-^ot\ pb{x,t,y) ) Jq 
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Proof. It clearly follows from ip G H'^'^ that V G L'^, hence V\ := VIa G L^- Therefore, 
by Fubini's theorem, J^^y Va{s) ds < oo for almost all x, y & W^, as desired. We claim 
that 



lim 



exp 



[-jlVixs)ds) 



1. 



exp(-£y(7f^)d.) 
Since one has, as a consequence of F o G (in this proof stands for L^(R; 

lim — / V^(7s) ds = — limt / V{x + {y — x)s) ds = 0, 
Jo Jo 

we just need to prove that 

exp(^- J V{xs)ds^ 



limE.^'^ 
t-*o ^ 



In fact, we observe that for any constant e > 0, the continuity of the paths of x. implies 
that 



lim 



qP^'^(|x. -7.^«|loo >£) =0, 



where L°° stands for L°°{[0,t];'R'^). Moreover 



exp ^ / —V{xs) ds^ — 1 



exp 



+ Ef'^ exp f -V{xs) ds^ - 1 



V{xs) ds) - l\\x. - < e 



\x. - 7.^^|l°° > £ 



Using the elementary inequality — 1| < 1 V e^, we can write 



exp J -V{xs) ds) - 1 < 1 V exp -V{xs) ds^ 



hence 



lim 

t^o 



< limEf'^ 

t^o ^ 



X- - 7.^^|lo° > £ 
X- - 7i^^|L'x. > e 



exp ^ ^ -V(xs) ds^ 

exp ^ y" -y(xs) ds) 
< lini (\ V e*(-'^^^)) - 7."^^^- > e) = 0. 

Similarly, using the elementary inequality — 1| < el^l — 1, we have 
exp ^ y ds^ - 



lim 



X- - 7i^^|L°° < e 



< limEj"''^ 



exp 



y(x,) + y(7,)|o!s 



X- - i^'^\L°^ < £ 



< lime^^e*l^°^^^li^i -1 = 0, 
~ t-^o 
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where we have used the following immediate consequence of the hypotheses: for any 
e > there exists 5s > such that 

rt 

\yils^) - yixs)\ ds = \V o -f'^y - V o x\li < 5e 

whenever \x. — ")f ^\L°° < £• This concludes the proof of the claim. Assertion @ now 
follows immediately from the claim just proved and the previous proposition. Assertion 
© follows by 



exp ( /g -V{xs) ds 



for t ^ 0, hence 



t^o t 
1 

t^o t 



Pb{x,t,y) 



lim — log E^'^ 



exp 



1 + oit) 



(V(y)-V(x))) 



-V{xs) ds 



lim — log ( exp 



-V{^s)ds) +o{t) 
lim - log exp (^t J —V{x + {y — x)s) ds^ 



-V{x + {y — x)s) ds, 



which finishes the proof. 



□ 



Remark 6. If x. is Brownian motion, Aizenman and Simon ^ proved that V is of Kato 
class if and only if 



lim sup E^ 



\Vixs)\ds 



0, 



without assuming that V is lower bounded. Using Kasminskii's lemma (see e.g. jl8j l it 
is then immediate to deduce that if V is of Kato class, then one also has 



(6) 



lim sup E^ 



exp 



\V{xs)\ds 



1. 



Moreover, Chung and Zhao proved that © continues to hold for any Hunt process x., 
under the only assumption that V is of Kato class (see Proposition 3.8). We also 
refer to [S] and [7] for related results on generalized Schrodinger operators and forms 
and associated Markov processes. 

4. Reconstruction of V 

In this section we assume that the hypotheses of Lemma|2land Proposition|3are satisfied. 

We shall then show that the transition probabilities of {xt,¥'^) determine the X-ray 
transform of Va. (in the sense of, e.g., ^1), which in turn yields Va. by a Fourier transform 
argument. In particular, from equations Q and © of Proposition 13 it immediately 
follows that 



F{x,y) := [ V{x + 
Jo 



{y — x)s)ds, x,y £ dA 



is determined by the transition probabilities of (xt,P^'). Moreover, one immediately 
recognizes that F is the X-ray transform of Va- Let us fix some notation: we represent 
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a line 7 as a pair 7 = (w, z), where u £ S"~^ is a unit vector in the direction of 7 and 
z G 7 n a;-*-. Then the hne integral f{x) dx is denoted by 

f{^)=f{oJ,z)=P^fiz). 

As argued before, since we know that V\ £ L^, Fubini's theorem implies that for each 
uj G S*^"^, PujVa{z) is defined for almost all z £ uj^. Moreover, we have, for p £ uj^, 

Va{p) = I PMz)e'^'''^ dz 



(a result often called slice-projection theorem, see e.g. ^l]). One uniquely recovers V\ 
by taking the inverse Fourier transform of Va(p)- Summarizing, we have proved the 
following 

Theorem 7. The restriction of V to the domain A can he uniquely reconstructed from 
the transition probabilities of (xj,P^). 

Remark 8. If our only aim were to reconstruct the function V\ from its X-ray or 
Radon transform, even more generality could be allowed, up to the situation where V is 
a distribution. For results on inverting the Radon transform of a distribution, see e.g. 

m, M- 

5. Reconstruction of the drift 

As in the previous section, we assume that the hypotheses of Lemma Inland Proposition 
Elare satisfied. 

Set u{x) = e'^^^'^~'^^y\ where y is any (fixed) point on the boundary of A. Then 
u satisfies the following Dirichlet boundary value problem for a second-order elliptic 
operator: 



(7) 



la'^Ux.Xj + b^Uxi = V{x)u, X £ A 



This is easily seen as a consequence of the definition of V and of the following simple 
calculations: 

Ux, = Ulpx^ 
Ux.xj = Uxjtpx, + Ulpx.Xj = u{lpxj'>Pxi + i^XiXj) 

where the last step is justified by 

a'^^x.i^x, = (aVV,VV) = {a^'^V^,a^'^V^) = la^/^VVp. 

If (O is uniquely solvable, then we are able to recover ip{x) for x G A uniquely. In fact 
we have: 

Proposition 9. Suppose a^^ are dijjerentiable and V £ L'^{A). Then there exists a 
unique solution u £ H'^'^{A) of the Dirichlet problem Q). 

Proof. Since ip £ C^{A), as follows by Sobolev embeddings, then f{x) := e^^^^~'^^y^ £ 
C^(A) C F^'^(A). Moreover, Lb,a is strictly elliptic and b = Vip £ C(A), hence b is 
bounded. We can now appeal to Theorem 8.9 of J3]) which yields the existence and 
uniquess of a solution to 0, as claimed. □ 
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Remark 10. Using more general results on elliptic PDEs, one can remove the unpleas- 
ant assumption of V being bounded, at the cost of added technicalities. In particular, 
using the existence and uniquess results of j2J , one can replace V G L°°{A) by 5 € L'^/^ 
in the hypotheses of Proposition |^ where 

The details (mostly calculations) are left to the reader. The assumption 1/ > is 
used to guarantee that the spectrum of the operator „ — V{x) (considered between 
appropriate spaces of integrable functions) does not contain zero. If we are willing to 
accept this level of generality, sacrificing a bit of concreteness, then we can dispense with 
the assumption of V being positive, and simply assume that zero is not an eigenvalue of 
I-'b,a ~ y{x). For further details we refer to jl3j and |U, where a Fredholm alternative 
for this type of operators is established. 

If a is the identity matrix, hence c = Vt/', we can obtain stronger results. In particular, 
the Dirichlet problem Q reduces to the Dirichlet problem for the time-independent 
Schrodinger operator with Hamiltonian — + V: 



(8) 




= V{x)u, X G A 



for which there exists a rich literature. We can apply, for instance. Theorem 4.7 of 
(see also ^). We shall denote by ta the first exit time of Brownian motion from the 
domain A. 

Theorem 11. Assume that V is of local Kato class, A is hounded and regular, and 
(9) x^W exp - ^ V{ws) ds 

is bounded. Then there exists a unique solution u G C(A) of the boundary value problem 

m- 

Proof. The conditions on V and A are needed in order to apply the above mentioned 
results of ^J. Moreover, since ip ^ C"^ a.s follows by Sobolev embeddings, and thus 
f{x) := e^^^)-'^^^) G C{dK), part (iv) of Theorem 4.7 in ibid, ensures that 

u{x) = [exp ( - j^^^ V{ws) ds)f{wr^) 

is the unique solution of (— + V)u = such that u G C(A) and u{x) = e^^^^~'^^y^ on 
dA. □ 

Remark 12. A simple sufficient condition guaranteeing that V is of Kato class is 
V G LP/2, p > d/2 (see e.g. p. 233). Therefore, if ^ G HP''^ with p > d, V is 
of Kato class. On the other hand, we were unable to find simple sufficient conditions 
ensuring that ® is bounded. Let us mention, however, that each of the following 
analytic conditions is sufficient: 

(i) Ttl dt is bounded; 

(ii) -iA + y > (in the sense of operators), or equivalently: 

(iii) The spectrum of — |A + F is contained in ]0, +oo[, 

where we have denoted by Tt the semigroup generated by —^A + V. For more informa- 
tions see jTj and p. 126. 
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Once u is obtained, one recovers ifj immediately, and hence c. We have proved the 
following result 

Theorem 13. Assume that the boundary value problem ^ admits a unique solution. 
Then the transition probabilities pi,+c{x,t,y) for x, y G h.^, t G [0,T], T > 0, determine 
c uniquely. 

Proof. It is just a combination of the previous steps. In particular, one proceeds as 
follows: 

(1) Obtain the X-ray transform of V\ from the transition probabilities ph^cix,t,y); 

(2) Invert the X-ray transform obtaining Va; 

(3) Solve the elliptic PDE ((TJ obtaining 7p{x) = logn(x) + ip{y)] 

(4) Obtain c = aVtp. 

□ 



6. Some extensions 

6.1. a^^c not a gradient field. It is clear that our approach strongly relies on the 
assumption that the a~^c is a gradient field. When this is not the case, we can only give 
a rather involved sufficient condition to reduce the problem to a more tractable one. We 
shall assume for simplicity that x. is a Lcj diffusion, without knowing a priori that c is 
a gradient field. We also assume that the technical assumptions introduced so far are 
in place when needed. 

Proposition 14. Let f : —)■ be a diffeomorphism mapping A into a bounded 
domain, and such that 

(10) [Vf{f-Hx))Vf{f-\x)r] -'L,jf{f-\x)) = VV(x) 

for some : M.'^ ^ R. Then the transition probabilities of (xj,P^) uniquely determine 
the drift c. 

Proof. Ito's formula for fi{xt), the i-th component of /, gives 

fi{xt) - fi{xo) = / {Vfi{xs),dws)+ / Lcjfi{xs)ds, 
Jo Jo 

that is 

f{xt) - f{xo) = / \7f{xs)dws+ / Lcjf{xs)ds, 
Jo Jo 

or equivalently, defining y. = f{x.) 



{ys)) dws- 



(11) yt = y+ f L,jf{f-\ys)) ds + f Vf{f-^ 

Jo Jo 

The hypotheses imply that one can reconstruct the drift and the transition probabilities 
of the diffusion ()11|). But this is enough to recover the transition probabilities of x. as 
well, as the following obvious identities show: 

P(xt = y|xo = x)= ¥{f-\yt) = y\f-\yo) = x) = ^{yt = f{y)\yo = f{x)). 

It is well known that the transition probabilities of an Lc,/-diffusion determine the drift 
coefficient c uniquely. □ 
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6.2. Unbounded A. We assumed A to be bounded in order to obtain existence and 
uniqueness of solutions for the Dirichlet problems ((TJ and 0. However, in some cases 
this assumption can be relaxed. For instance, imposing enough boundedness of the 
coefficients in ((JJ, one can obtain existence and uniqueness results without assuming 
that A is bounded. 

Proposition 15. Assume that a*-' are differentiable, e~'^' \\a{x)\\ — > for \x\ — > oo, and 
V is non-negative (in the generalized sense). Then there exists a unique solution of Q). 

Proof Let v{x) = n(x)e^(^) = e'^(^)+^(^), ip{x) = -\x\^. Since G H"^'^, hence £ 
by Sobolev embedding, ip is bounded on A, and so is u. Moreover, it is immediate to 
show that V € H'^'^{A). Then one has 

and Vw = (Vtp + Vip)v, 

Lb,aV = {Lb,a{i^ + ip) + \a^^^V{'^ + ip)f)v 

= tr(aV'xx) + {h + 2aV<f, W) + |a^/^VV'P + Lb,af + \a^^^V^\'^y, 

therefore 

Lb-2aW<p,aV = Lb,aV " {^aVif, Vv) 

= Lb^aV — {2a'V(p, Vip + Vip)v 

= Lb,aV - {2aVLp,V7p)v -2\a^/^Vip\'^v 

= [I tr (aV'x..) + {b, V^) + la^/^VVp + Lb^aV - W'^V^\^)v 
= {y{x) + d{x))v, 
where d := Lb^af — \a^^'^'Vip\'^ . That is, v solves the equation 

(12) Lb-2aVv,aV = {y{x) + d{x))v, 

with boundary condition v{x) = /(2;)e'^(^) on dh. One can apply now Theorem 8.9 of 
jl3j to determine existence and uniqueness of a solution of (|12p in H^'"^. In fact, if L^ ^ 
is strictly elliptic, so is also Lb-2a'Vtf,a, and the continuity and growth condition on a 
imply that b — 2aS/ip is bounded. Finally, there is no loss of generality in assuming that 
y > implies V{x) + d{x) > 0: if it were not so, we could make \d\ arbitrarily small by 
using as cut-off function (p{x) = e~'^l^l , without altering the previous results. Then the 
unique solution to is given by u{x) = t;(x)e~'^(^'). □ 

In the case of unit diffusion a stronger statement can be made, as follows by the 
results in Chapter 5 of jllj . 

Proposition 16. Let d > 3, V £ L^i-^) and Kato class. If uj ^ oo in A, then the 
solution of 10) is given by 

it/(x)=E^ exp(^-^ V{ws)ds^f{wr^) 

Proof. As before, ip € H"^'"^ implies ip £ C^, hence that ip is bounded on dA, and also 
that fix) := e'/'(^)-'^(?') G L^{dA). Then by Theorem 5.18 and 5.19 of [HI we obtain 
that Uf solves (— + V)u = and Uf € Cf,(A). Finally, Uf satisfies the boundary 
condition as an immediate consequence of its definition and continuity in the closure of 
A. □ 
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6.3. One-dimensional case. Using the considerations of the previous subsections it is 
possible to give a solution, at least in some special cases, to the problem posed in [2] of 
reconstructing the drift of a one-dimensional diffusion with a = 1. In particular, let us 
assume that the transition probability density pi{x,t,y) of the diffusion 

dxi{t) = ci(xi(t)) dt + d'Wi{t) 

is known for x, y €]0, t € [0, T], with T a fixed positive constant. As before, our aim 
is to determine ci(x) for x G [0, 1]. 

Define an M^-valued diffusion as weak solution of the following SDE: 

f dxi (t) = ci{xi{t)) dt + dwi (t) 
(13) <^ 

[ dx2{t) = C2(x2{t)) dt + dW2{t), 

where C2 is a smooth, bounded function, and wi, W2 are standard independent Brownian 
motions. In more compact notation, we can write 

dx{t) = c{x{t))dt + dw{t), 

with c{xi,X2) = (ci(xi), 02(^2)), w = {wi,W2)- It is clear that one can recover the 
transition probabilities of xi from those of x, since xi and X2 are independent. 
Trying to reconstruct the vector field c (which is trivially a gradient field: c = V^/^, 
■0(xi, X2) = Jq^ ci(s) ds + J^^ C2{s) ds), one is faced with two problems: A =]0, l[xM is 
unbounded, and Va is not in L^. One can try, however, to consider the problem of drift 
reconstruction for y. = f{x.), with / a diffeomorphism. 

Proposition 17. Assume that there exist C2 : M ^ M and a diffeomorphism f :M? 
satisfying Iil0\) and such that the corresponding Dirichlet problem admits a unique 
solution on /(A). Then ci can be uniquely reconstructed from the transition probabilities 
of xi{-). 

Proof. The transition probabilities of xi(-) uniquely determine the transition probabil- 
ities of x{-) outside the (unbounded) domain A =]0, l[xM. If / exists such that (|1U|) 
holds, then the the problem of reconstructing the drift of y(-) := f{x{-)) is well posed, 
and it is equivalent (under the technical assumptions introduced in the previous sections) 
to the solvability of the Dirichlet problem ((T)) on the domain /(A). Therefore, assuming 
the latter problem admits a unique solution, this yields the transition probabilities of 
y(-), hence those of x{-) because / is a bijection. As already observed, the transition 
probabilities of a diffusion with generator Lcj uniquely determine c. □ 

It is clear that the last proposition is not constructive and simply gives sufficient 
conditions for the solvability of the problem of drift reconstruction in dimension one. 
As in the case of higher dimensional diffusions with a~^c not being a gradient field, 
these sufficient conditions seem difficult to check. However, since we essentially rely on 
the above described representation with the drift being a gradient field, this seems to 
be the best we can achieve by our present method. 
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